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Abstract. After summarizing the most interesting results in the calculation of synchrotron 
radiation in the Myers-Pospelov effective model for Lorentz invariance violating (LIV) 
electrodynamics, we present a general unified way of describing the radiation regime of LIV 
electrodynamics which include the following three different models : Gambini-PuUin, Ellis et 
al. and Myers-Pospelov. Such unification reduces to the standard approach of radiation in a 
dispersive and absortive (in general) medium with a given index of refraction. The formulation is 
presented up to second order in the LIV parameter and it is explicitly applied to the synchrotron 
radiation case. 



1. Introduction 

The appearance of detectable low energy effects arising from possible quantum gravity 
corrections to particle dynamics ^ has been recently the subject of intense scrutiny 
The most direct interpretation of such effects, though not the only one is related to the 
breaking of observer Lorentz covariance. In this way this subject becomes directly tied to the 
long time honored investigations, both theoretical and experimental, based upon the Standard 
Model Extension [3j concerning Lorentz and CPT violations [H]. Heuristic derivations of such 
effects from loop quantum gravity have been presented [ZllHl) which make clear that a better 
understanding of the corresponding semiclassical limit is required [HI- String theory has also 
provided models for explaining such quantum gravity induced corrections JHl- Also, effective 
field theory models to describe such effects have been constructed [TT] . 

The observation of 100 MeV synchrotron radiation from the Crab Nebula has recently been 
used to impose very stringent limits upon the parameters describing a modified photon dynamics, 
embodied in Maxwell equations that get correction terms which are linear in the Planck length 
(Hj. Such bounds are based on a set of very reasonable assumptions on how some of the 
standard results of synchrotron radiation extend to the Lorentz non-invariant situation. This 
certainly implies some dynamical assumptions, which have been recently tested by two of us 
(RM and LFU) J2] in the framework of the classical version of the Myers-Pospelov (MP) 
effective theory, which parameterizes LIV using dimension five operators, as described in the 
first paper of Ref. Ref. J2j summarizes a complete calculation of synchrotron radiation in 
the context of this model. This constitutes an interesting problem on its own whose resolution 
will subsequently allow the use of additional observational information to put bounds upon 



the correction parameters. For example we have in mind the polarization measurements from 
cosmological sources. The case of gamma ray bursts has recently become increasingly relevant 
|13j . although it is still at a controversial stage |14j. 

This paper is organized as follows. In Section 2, we provide some highlights of the results 
obtained in Ref. jl2j . and from there onwards we develop a formalism for a Lorentz violating 
electrodynamics based on extended constitutive relations for the D and H. In Section 3 we 
state the general structure of this theory, in Section 4 we introduce a parameterization for the 
constitutive relations, and in Sections 5 and 6 we summarize the main expressions for the fields, 
energy spectrum, power spectrum and polarization in terms of effective refraction indices. We 
close with a short summary in section 7. 



2. Synchrotron radiation in the Myers-Pospelov model 

This model parameterizes LIV using dimension five operators both in the matter and 
electromagnetic sectors. There is also a preferred frame four velocity V^, which is not a 
dynamical field. We choose to work in the rest frame = (1,0). As usual the model possesses 
observer Lorentz covariance, which means that the fields and the four-velocity transform as 
tensors when going from one observer frame to another. On the other hand, in each frame we 
violate particle Lorentz transformations; that is to say we have a collection of non-dynamical 
physical fields in the action, in analogy to the physics going on in the presence of an external 
magnetic field which violates particle rotation invariance, for example. 
In the rest frame, the modified Maxwell equations are (c = 1) 

dB 

V-B = 0, VxE+— = 0, 

ot 

dE 

V-E = 47rp, -— + VxB + e9o(-VxE + aoB))=4^j, (1) 

where the LIV parameter ^ is usually written as ^ = (,/M with ^ being a dimensionless number 
and M being a mass scale characterizing the Lorentz symmetry breaking, which is usually, but 
not necessarily, identified with the Plank mass. For the particular case of synchrotron radiation 
the modified dynamics for a particle of charge q in an external constant magnetic field in the 
z-direction produces the standard helicoidal orbits with a modified Larmor frequency loq 

r = ^(l-^vE+ ^fE^^ (v X B) ^ (v x B) , (2) 

depending upon another LIV parameter fj = —r]/M and with E being the energy of the particle. 
The minus sign is chosen to make contact with the notation of Ref.^. We restrict ourselves 
to circular orbits in the x — y plane with Larmor radius R = P/ujq and we use the standard 

1/9 

notation /? = |v| and 7 = (l — . When rewriting (3 in terms of the particle energy we 

obtain the modified expression 



(3) 



Due to de admixture of vector and axial vectors in the last Eq. the problem presents 
birefringence and the two independent modes of propagation correspond to those of definite 
circular polarization A = it, with refraction indices n\{^z) 



n\{z) = VT + z^ + Az, z = ^uj. 



(4) 



Such modes are identified via the corresponding Green functions, which are most easily calculated 
in the radiation gauge. It is interesting to observe that Eqs. (^Q) provide the exact expression 



47rS = E X B - X 5E/5t, 
for the Poynting vector. Also, the magnetic and electric field satisfy the relation 

B(a;, fi) = \/l + z^h x E(a;, n) — iz E(u;, n). 



(5) 



(6) 



in the radiation approximation, with n = r/r being the direction of observation. These two 
effects conspire to produce a positive definite expression for the power flux yielding the following 
expression for the average angular distribution of the radiated power spectrum 



'£P{T)\ 




dR 



m, A 



dn 



UJr, 



muJo, 



4^ ^/TT^ 



[\l3nx{Zm)JL{Wxm) + cote Jm{Wxm)\ 



(7) 
(8) 



Here, the average () is taken with respect to the macroscopic time T and J^, Jm denote the 
Bessel functions and their derivatives respectively. The argument of the Bessel functions is 
W^Am = fnn\{zm)l3sm9. We also have calculated the total averaged and integrated power 
radiated into the m-th harmonic 



m A=± 



J2^(2m nx{zm)f3) 



2p^nlizm) Jo 



dx J2m{x) 



(9) 

which clearly indicates the contribution of each polarization Pxm- The above result is exact in 
Zm and the parity-violating contribution has vanished after the angular integration. 

In the case of synchrotron radiation from Crab nebula, as well as from other objects like Mkn 
501 and GRB021206, one can estimate that the m corresponding to the observed radiation are 
very large, let us say in the range 10^^ — 10^^, with the respective values of m/7 varying within 
10"^'^ — 10^''. These properties motivate the large m and m/7 expansion of the Bessel functions 
involved in the spectrum. In a similar way to the standard case [TF, a first consequence of this 
approximation is the appearance of the cutoff frequency cjac = ^0 with 



l-(3\E)nl{zm) 



-3/2 



(10) 



This means that for m > rhxc the total power decreases as P\m ^ cxp(—m/rh\c). Within the 
same large-m approximation, the integrated power in the m-th harmonic can be expanded to 
second order in ^ yielding 
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where m^ 



Here K, 



p/q 



= 37^/2 and k(x) = x dy K^i^[y) is the so called bremsstrahlung function, 
denote the Bessel functions of fractional order. Let us notice the appearance of 



the combination ^u)m/^ = ^{Lo/Mp){m/^) as the expansion parameter in This is not 

necessarily a small number, which signals the possibility that such corrections might be relevant 
in setting bounds upon ^. This rather unexpected effect is due to the amplifying factor m/7. 



Similar results have been obtained in calculations of the synchrotron radiation spectra in the 
context of non-commutative electrodvnamics jl6j . 

Another possibility for observable effects due to is to look at the averaged degree of circular 
polarization 

_ (P+H-P-M) .... 

where P\{t^) is the total power distribution per unit frequency and polarization A, so that 
Px{uj) = Pmx/ojQ. The average here is calculated with respect to an energy distribution of 
the relativistic electrons, which we take to be N[E)dE = CE~^dE, in some energy range 
El < E < E2, chosen as E'l = and E2 ^ 00 for simplicity. The result is 

^e = ^^(^) n(p), (13) 



qB 



(p-3)(3p-i) (p + i)r(| + i)r(| + 



"^^^^ 3(3.-7) (^ r(| + M)r(| + ii) ' ^''^ 

Again, we have the presence of an amplifying factor in Eq. (|13() . given by {fiu/qB), which is 
independent of the form of Il{p) and not necessarily a small number. An estimation of this 
factor in the zeroth-order approximation {S, = = fj), which is appropriate in (|13j) . yields 
(fiuj/qB) = oj/{ujQ^) = m/7. The expression (|13)) is analogous to the well-known average of the 
degree of linear polarization Ulin = (p + l)/(p + 7/3), under the same energy distribution for 
the electrons. 

3. General structure of Lorentz violating electrodynamics 

Three paradigmatic examples of Lorentz violating electrodynamics are given by the effective 
theories proposed by Gambini and Pullin (GP) [7j, Ellis et al.fEMN'l|lUj. and Myers and 
Pospelov (MP) 11]. They can be written in the general form of Maxwell equations 

V • D = 47rp, V • B = 0, (15) 

V X E = - — , V X H = — + 4vrj, (16) 

with corresponding constitutive relations 

D = D(E,B), H = H(E,B), (17) 

which we next write in detail for each case, after reviewing the corresponding equations. Let 
us recall that the above equations (|15j) and (fT6|) imply charge conservation dp/dt + V • j = 0, 
independently of the constitutive equations (fT7|) . In an abuse of notation have denoted by ^ the 
electromagnetic LIV parameter for all models in the sequel. 

3.1. Gambini-Pullin Electrodynamics 
The Maxwell equations for this case are 

V-B = 0, V X (S + 2iV X S) + — = 0, (18) 



dt 

V-^ = Anp, V X {b + 2lV X b) - — = 47rj, (19) 



where the electric and magnetic fields are identified from the homogeneous equation as 

E = f + 2|V X <5, B = B. (20) 

From the inhomogeneous equations we obtain 

D = H = e + 2|V X (21) 

which together with the constitutive relations 

D + 2|V xD = E, H = B + 2|V x B, (22) 

leave the equations in the required form. These equations define the corresponding functions 
stated in (|17() . In momentum space 

D = K (e - 2i|k X E + 4|2 (k • E) k) , H = B + 2i^k x B. (23) 

1 + 4^2/j,2 V / 

The admixture of vectors and axial vectors in the constitutive relations precludes the parity 
violation exhibited by the model, together with the presence of birefringence. 

3.2. Ellis et. al. Electrodynamics 

In this case the modified Maxwell equations are 

(?B 

V-B = 0, VxE + — = 0, (24) 

5E 

V-E + U- — =4^Pe// =47r(p-u-j), (25) 
5E 9B 

VxB - (l-u2j — + UX — + (u-V)E = 47rje// = 47r(j + u(p-u-j)). (26) 

which, following the approach of Ref. ^0] and assuming that in momentum space u = /(a;)k, 
can be written in the form (|15lll6j) via constitutive relations which read 

H = B-/(u;)kxE, D = (l - /^(u;)^^) E + /2(u;)k (k • E) - /(a;)k x B. (27) 

Taking u as a vector, this model conserves parity and shows no birefringence. 

3.3. Myers-Pospelov Electrodynamics 

This case corresponds to the equations (^. From the last one we can infer the constitutive 
relations 

H = B - idoB, D = E-|5oB, (28) 

which produce 

V • E =V • D, (29) 

leaving the third Eq. in desired form. Similarly to the GP case, this model violates parity. 
In momentum space Eqs. ()28() become 

H = B + iiojB, D = E + i^LoB, (30) 



The above constitutive relations in the three representative models involve linear relations 
among the fields and can be summarized in momentum space as local relations 

Di{uj,k) = aij{LO,k)Ei{uj,k) + pij{Lj,k)Bj{uj,k), 

Hi{uj,k) = Pij{u;,k)Bi{u;,k) + aij{uJ,k)Ej{uJ,k), (31) 

where the corresponding momentum dependent coefficients can be read from the equations H23|) . 
()27|). and (|3U|) . Equations pij) are the most general expressions in which any pair of linear 
constitutive relations can be ultimately written, because one must be able to solve the fields 
D, H in terms of E, B. 



4. Parameterization of the constitutive relations 

Let us consider Maxwell equations in momentum space 

k-B(tJ,k)=0, kx E(a;,k) =a;B(w,k), (32) 

ik • D (u;,k) = 47rp(u;,k) , zk x H (w, k) = -icjD (w, k) + 47rj (a;, k) . (33) 

Here we discuss the vacuum situation where the non trivial constitutive relations arise because 
of LIV effects. Let us take into account corrections up to second order in the LIV parameter 
and let us assume that we are in a Lorentz frame where we demand invariance under rotations. 
This would correspond to the rest frame = (1,0) in the MP model, for example. We can 
always go to an arbitrary frame by means of an observer Lorentz transformation. In this way 
we have the general expressions 

ttij = ao6ij + iai^eirjkr + a2i^kikj, pij = poSij + ipiieirjK + P2£.^kikj, (34) 
Pij = PoSij + iPiicirjkr + P2i'^kikj, Gij = a^Sij + iaiicirjkr + CF2i'^kikj, (35) 

where aA, Pa-, PAi^^A-, = 0, 1,2, are escalar functions depending only upon oj, k = |k|, and ^. 
The property k • B = sets /?2 = /O2 = effectively. In vector notation we then have 

D = (ao + 02^:^1^) E + (po + iaiuj^) B + i (pi - ia2w|) | k x B, (36) 
H = (cto + a2k'^i'^) E + (/3o + i(Jiioi) B + i - ia2U)i) | k x B, (37) 

where we have used the second Eq. (|32j) together with (k • E) k = w (k x B) + fc^E. Next we 
substitute in Eqs. (|33|l to obtain the corresponding equations for E and B. The result is 

i (ao + oi2k^i^^ (k • E) = 47r/9, (38) 



i {ctQ + a2k'^i'^^ wE + i I3q + i {cri + pi) + 02!^ o;^ 



k X B 



(ao + po)uj + i (aiuj^ - Pik^) I] B = 47rj (w, k) . (39) 



Let us rewrite the inhomogeneous equations in the compact form 

iP{k-B) = Attp, (40) 
iujPB + iQkxB + RB = 47rJ(a;,k), (41) 

by defining 

P = ao + a2f^A;2, Q = pQ + i (ai + pi) + a2S,'^uj^ , R= (^Pik"^ - aiuj"^^ ^ + i {ao + po) uj . (42) 
Now we have only three independent functions which depend on oj and k. 



5. Potentials and fields 

Starting from (p^ . (I^U]) and (|1T|) we introduce the standard potencials ^, A 

B = ik X A, E = iioA-ik^. (43) 

Substituting in the inhomogeneous equations we have 

-wk-A + A;2$ = 4iTp/P, (44) 
(^k'^Q-uj^P^A + iR{kxA)-[Q{k-A)-Puj^]k = Anj. (45) 

We can easily fix two convenient gauges as follows. 

5.1. Generalized Lorentz gauge 

In the standard situation the Lorentz gauge in momentum space is k - A = to^. Here we choose 

k-A = {P/Q)uj^, (46) 
in such a way that Eqs. (|45j) and (|44|1 reduce to 

(Qk^ - Pu^^ A + iR (k X A) = 47rj (47) 
(Qk^ - Puj"^) ^ = An {Q /P) p, (48) 

respectively. 

5.2. The radiation gauge 

As usual we set k • A =0 and we obtain the equations 



$ = 47r(^rPj p, (49) 
Qk^ - Puj^) A + iR {kx A) = 47rj - 4TTUjpk~^k. (50) 



Eq. (|49|1 exhibits the scalar potential as a static contribution, which subsequently does not 
contribute to the radiation field. The use of current conservation in (|5nj) allows the introduction 
of the transversal current, consistently with the chosen gauge. That is to say we have 

(Qfc2 - Plj^^ a + ikR (k X a) = Att [j- (j Att^t- (51) 

The presence of birefringence depends on the term proportional to kR, where we have a crucial 
factor of i. This makes clear that a possible diagonalization can be obtained by using a complex 
basis, which is precisely the circular polarization basis. In fact, decomposing the vector potential 
and the current in such basis 

A = A+ + A-, jT=jT+j?, (52) 
and recalling the basic properties 

kxA+ = -iA+, kxA" = iA-, (53) 
we can separate (|?T|) into the uncoupled equations 



Qk^ - Poj^ + XkR A^ = A-Ki^, A = ±1. (54) 



6. Green functions in the radiation gauge 

To our aims it is convenient to work in this gauge, where the polarized Green function is defined 
by 

(Qk'^-Puj^ + XkR)G^iLj,k) = l. (55) 



For each polarization, the above leads to a cubic equation in k, which determines the dispersion 
relations for the propagating photon. In fact one of the poles is due to the Pi factor in the 
expression for R in Eq. (|42j) . This pole is at A; and therefore its contribution can 

be neglected in our effective theory valid for k <^ We defer for future work a detailed 

characterization of the solutions and proceed here on the basis of general properties of light 
propagation, which constitute a generalization of the specific models already considered. We 
characterize the propagation mode corresponding to each polarization A by a refraction index 
n\{(jj) to be read from the appropriate dispersion relation in such a way that n\{uj) = k\{io)/io. 
The general form of the polarized Green function is 

G\uj,r) = -^F^(u;)e*-"^(^)^ (56) 
47rr 

where n = r/r is the direction of observation. The generic form for the total Green function in 
the far-field approximation is 

Glit {u, r) = i [(5^'=-n^n'= + ie'^^h'^ G+ (w, r) + (^6'''-n'n'' - ie^'^r) G" {u, r)] . (57) 

Note that from the birefringent case we can go to the non-birefringent one by taking n-^{uj) = 
n_(u;) = n{Lo), in which case F~^{uj) = F^{uj) = F(uj) and the Green function H57|) becomes 

Gilt i^, r) = 7^ f 5^^-n*n^) F (uj) e*""(")'-. (58) 
47rr V / 

In the following we analyze the main consequences of an electrodynamics characterized by a 
Green function of the type (|57|) . To warrant that the fields are real, it must be 

[G+(u;,r)]* = G-(-u;,r). (59) 

This implies the relations 

nXiio) = n.{-uj), [F+{uj)]* = F-{-uj). (60) 

For a birefringent medium the real and imaginary parts of the refraction index for circular 
polarization components can contain both w-even and w-odd terms, provided that they satisfy 
Eq. ()6U() . In the case of a non-birefringent medium the real part of the refraction index is even in 
uj, and the imaginary part is odd. We can see that the refraction indices for the Myers-Pospelov 
theory, Eq. @, satisfy these requirements. 

The Green functions must also be consistent with causality, which means that they cannot 
have poles in the uj upper half plane. This leads to generalized Kramers-Kronig relations, which 
for a dispersive refraction index imply the existence of an imaginary part. This means that 
these media are necessarily absortive. If we assume that these theories are valid in a range of 
frequencies where the effective medium can be considered as transparent, then in the case of 
birefringent theories we have a real refraction index with both a;-even and w-odd terms, while 
for non-birefringent theories we have only w-even terms. 

Next we turn to the calculation of the angular distribution of the power spectrum. The full 
vector potential is given by the superposition 

A{u;, r) = -J2 F\oj)e'^^^'3\uj, k,) = A+ + A_. (61) 

A=±l 



where = nx{uj) Loh and n = r/r. From here we can compute the electric and magnetic fields 
B(u;,r) = V X A(w,r) = - ^ AnAFV"^'^''j^(w,kA) = u;(n+A+ -n_A_), (62) 



A=±l 



E(u;, r) = iioA{io, r) = — ^ F^e*"^'^l^(c^, kA) = iio (A+ + A, 



A=±l 



The Poynting vector 



(63) 



S(a;,r) = — E(-u;, r) x H((j, r) 
4tt 



7~ X! (^0 + ^'^1'^'^) '^A + «Acro [A+ (Aw, r) x A_ {-Xuj, r)] 



A=±l 



+XPiiujnl [A+ {Xuj, r) • A__ {-Xuj, r)] n| , 
allows us to compute the angular distribution of the energy spectrum 



(64) 



dilduj 



27r 



n- [S {uj,r) + S{-uJ,r)], 



'Z{uj) A_ (-(J, r) • A+ {uj, r) + Z{-uj) A {uj, r) • A+ {-uj, r)] , (65) 



where 



Z{uj) = ^ J ifjo (Aw) + 

A=±l 



Pq{Xuj) + iA 0"! (Aw) CJ,^ 



nx{Xuj) + f3i{Xuj)iujnl{Xuj)\ . (66) 



Using Eqs. (|6()|) and (|61|). together with the general property [j'*'(w,k)]* = j '^(— w, — k*), which 
allows for complex refraction indexes, we finally get 



4^ = ^12 ZiXuj) F\uj) 
dClduj 87r2 ^ ^ ^ ^ ^ 



A=± 



-2(^Im[n;\(<^]r 



j^(w,kA(w)) •j\w,kA(w)), (67) 



where the vectors kA(w) were defined after Eq. (|61j) . As expected, absorption is present through 
the non-zero imaginary part of the refraction indexes. Eq. (|67|l is related with the angular 
distribution of the radiated power spectrum by 



dQduj 

where T is a macroscopic time, which leads to 



dT 



d^P{T) 
dwdVL 



(68) 



d^P{T) 
dClduj 



UJ 

8^ 



E 

A=± 

oo 



Z(Aw) 



FHuj] 



-2ujlm[nx{u!]r 



/oo 
dT e-*-^j^(r + T/2,kA(w)) Pt jV(r + r/2,kA(w)) 
-oo 



(69) 



where P^^ = 1/2 ( 6kr — n^nr + iXekirfii) are the circular polarization projectors. 

In the sequel we consider synchrotron radiation in the situation where the vacuum absorption 
can be neglected so that the refraction indices are real and Eq. (|6U() reduces to n^{uj) = n~{—uj). 



Also we assume that F (uj) and Z (lo) are real functions. Using the current corresponding to the 
circular orbits in we can compute the cycle average over the macroscopic time T obtaining 



'd2p(r)\ 



/ A=±l m=0 



dPr 



dn 



(70) 



with 



— -Z[\oj) 



dn 



Svr 



F\uj) 



[X(3nx{u)J:^{W^x) + cot9 MWrr,x)]\ 



where W^x = mn^{Lo) (3 sin 9. The total averaged power radiated in the m^^ harmonic 



(71) 



^2 /o2 



2R 



A=±l 



J2^{2m(3nx) 



1 - P^nl f-impr^x 



2(5^nl Jo 



dx J2m{x) 



(72) 



where the contribution of each polarization is exhibited. 

In the case of the Myers and Pospelov effective theory the non-null coefficients in the 
constitutive relations (|H(:ii;-i7jl are 

ao = /3o = 1, cro = Po = ii^, (73) 



while the remaining functions are given by 



F\u^) = ^^^=, Z(^) = 2Vl + (e^)2. (74) 
From the general expressions obtained in this section we recover the results of Section 2. 



7. Final remarks 

We have presented a compact review of the most interesting results for synchrotron radiation 
in the Myers-Pospelov LIV electrodynamics. This is one of several existing models of LIV 
electrodynamics that also include those proposed by Gambini-Pullin, and Ellis et. al., for 
example. Motivated by such diversity we have developed a general approach based on generalized 
constitutive relations, set out in terms of the standard methods for describing radiation in a 
dispersive and absortive medium with a given index of refraction. The resulting formalism has 
been explicitly applied to synchrotron radiation in the case of transparent media. 
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